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ABSTRACT 

We construct an analytic formalism for the mass function of cold dark matter 
halos, assuming that there is a break in the hierarchical merging process. Ac- 
cording to this broken-hierarchy scenario, due to the inherent nature of the gravi- 
tational tidal field the formation of massive pancakes precedes that of dark halos 
of low-mass. In the framework of the Zel'dovich approximation which generically 
predicts the presence of pancakes, we first derive analytically the conditional 
probability that a low-mass halo observed at present epoch was embedded in 
an isolated pancake at some earlier epoch. Then, we follow the standard Press- 
Schechter approach to count analytically the number density of low-mass halos 
that formed through anti-hierarchical fragmentation of the massive pancakes. 
Our mass function is well approximated by a power-law dN/dM = M~ l in the 
mass range ltfh~ l M Q <M< 10 w h' 1 M Q with the slope I = 1.86 shallower than 
that of the currently popular Sheth-Tormen mass function I = 2.1. It is expected 
that our mass function will provide a useful analytic tool for investigating the 
effect of broken hierarchy on the structure formation. 

Subject headings: cosmology: theory — large-scale structure of universe 



1. INTRODUCTION 

In the cold dark matter (CDM) paradigm the gravitationally bound objects made of 
dark matter particles (dark halos) are supposed to form through hierarchical merging process. 
Press & Schechter (1974, hereafter PS) devised for the first time an analytic formalism to 
evaluate the mass distribution function of CDM halos that formed hierarchically. Later Sheth 
& Tormen (1999, hereafter ST) refined the PS mass function by taking into account the non- 
spherical collapse. The ST mass function has been tested against many N-body simulations, 
showing good agreements (e.g., Reed et al. 2003). This success of the ST mass function 
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implies the validity of the hierarchical merging scenario since it was originally derived under 
the assumption that the CDM halos form through hierarchical merging. 

Yet it is premature to assert that the formation of CDM halos is always hierarchical 
over the entire mass range given the fact that the validity of the ST mass function was 
rather limited to the relatively high- mass section (> l^h^M®). Results from N-body 
simulations still suffer from the lack of information on the low-mass section (< lO lo /i- 1 M ) 
due to the resolution limit. 

Very recently, Mo et al. (2005) came up with a new halo-formation scenario where 
the formation of low-mass halos is somewhat anti-hierarchical, preceded by that of massive 
pancakes. According to their model, the preheated medium caused by the gravitational 
pancaking effect suppresses the star-formation sufficiently in the low-mass halos, which can 
explain the observed low Hi-mass as well as the faint-end slope of the galaxy luminosity 
function. 

Although the work of Mo et al. (2005) was focused on explaining the suppression of 
star-formation in the low-mass halos, we note that if the formation of low-mass halos was 
indeed preceded by the formation of massive pancakes, what is suppressed should not be 
only the star-formation but also the formation of low-mass halos itself. To take into account 
the break in the hierarchical process and to predict the abundance of low-mass halos more 
accurately, it is desirable to have an analytic model for it derived from first principles. 

Our goal here is to construct an analytic model for the low-mass halos that form through 
anti-hierarchical fragmentation of massive pancakes. To achieve this goal, we adopt the 
Zel'dovich approximation as a simplest footstep which generically predicts the formation of 
pancakes, and we follow the standard PS approach to count the abundance of dark halos 
as a function of mass. The hypotheses, the mathematical layout and the predictions of our 
model are presented in §2, and the summary and discussion of the final results are provided 
in §3. 

2. ANALYTIC FORMALISM 

In the Press-Schecther formalism (Press & Schechter 1974, hereafter, PS), an isolated 
dark halo (a halo just collapsed) of mass M forms from the regions in the density field whose 
average density contrast 5 = Ap/p (p: the mean mass density of the universe) smoothed 
on the mass scale of M reaches some critical value, 5 C . The value of the critical density S c 
is approximately 1.68, which depends very weakly on the background cosmology and the 
redshift (Kitayama & Suto 1996). 
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The Gaussian probability distribution of the linear density contrast smoothed with the 
sharp k-space filter on the mass scale M is given as 

"In k c 
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a 2 (M) = j C A 2 (k)d\nk, M = 6ir 2 pk~ 3 , (1) 
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where a(M) is the rms density fluctuation, and A 2 (k) is the dimensionless power spectrum. 
Throughout this Letter, we use the power spectrum of the concordance ACDM cosmology 
with Q A = 0.7, Q m = 0.3, Q b = 0.044, h = 0.7 (Bardeen et al. 1986). 

In the Zel'dovich approximation (Zeldovich 1970, hereafter, ZEL), the mass density is 
given as 

p= (i_A 1 )(l-A 2 )(l-A 1 )' (2) 
where Ai, A2, A3 (with Ai > A 2 > A3) are the three eigenvalues of the deformation tensor dij 
defined as the second derivative of the perturbation potential = didjQ. Doroshkevich 

(1970) derived the joint distribution of the three eigenvalues: 
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P(Ai, A 2 , A 3 ) = exp (--± + ^-|J (Ai - A 2 )(A 2 - A 3 )(A X - A 3 ), (3) 

with Ii = Ai + A 2 + A 3 and I 2 = A X A 2 + A 2 A 3 + A 2 A 3 . 

Equation (2) implies that the mass-density diverges along the direction of the major 
principal axis of the deformation tensor if the largest eigenvalue reaches unity and the inter- 
mediate and the smallest eigenvalues are less than zero. In other words, an isolated pancake 
(a two-dimensional object just collapsed only along the first principal axis) of mass M forms 
if the following condition is satisfied: Ai = A c , A 2 < 0, A 3 < with A c = 1 on the mass scale 
M. 

Now that the conditions for the formation of isolated halos and pancakes are specified, 
we would like to find the probability that a halo at present epoch was embedded in a 
pancake at some earlier epoch before it formed. For this, it is required to have the joint 
distribution of the linear density and the three eigenvalues of the deformation tensor on two 
different scales on two different epochs. Let 5 be defined at present epoch on the galactic 
mass scale M g , and let Ai, A 2 , A3 be defined at some earlier epoch of redshift z > on the 
larger mass scale M p > M g . The rms density fluctuation at redshift z on mass scale M p is 
given as <t(M p , z) = b(z)cr(M p ) where b(z) is the growth factor of the linear density that is 
normalized to satisfy 6(0) = 1. Since the growth factor as well as the rms density fluctuation 
is a decreasing function of z, we have a(M p , z) < cr(M g ). From here on, we use the notations 
of a'^, A^, A2, A3 to represent the rms density fluctuation, the deformation tensor and its 
three eigenvalues at redshift z on the mass scale M p . 
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To derive the joint distribution of 5 and A'^A^A'g, we first derive the multivariate 
Gaussian distribution of 5 and the 6 independent components of the symmetric tensor, 
d'^. Rotating the frame into the principal axis of d'^ and using the fact that 5 is invariant 
under the axis-rotation, we derive analytically the joint distribution of 5 and A' 1? A' 2 , A' 3 



p(S, Ai, A' 2 , A' 3 ) 
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where a\ = a 2 - a' 2 , I[ = X[ + A 2 + A 3 , and I[ = A;A 2 + A' 2 A 3 + X[X' 3 . 



Through equations and (1) and (4) we find the probability that a halo of mass M g 
observed at present epoch was embedded in an isolated pancake of larger mass M p at redshift 
z with the help of the Bayes theorem: 



p (X[ = A C ,A 2 < 0,A 3 < 0|5 > 6 C ) 



p(5 > S c , X[ = A c , A' 2 < 0, A 3 < 0) 
P{6 > S c ) ' 
/^/_° 00 dA^ 00 d.5p( ( 5, A', = A c , A' 2 , A 3 ) 
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In equation (5), the integration over 5 can be readily evaluated 
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This probability (eq.[5]) will allow us to determine the most-likely epoch when the formation 
of pancakes precedes that of low-mass halos, and the typical mass scale for the formation of 
pancakes as well. 

For comparison, we consider the probability that a halo of mass M g observed at present 
epoch just formed hierarchically at redshift z which is approximately given as (Bower 1991; 
Lacey & Cole 1994): 



p(S" = 5 C \5 > S c ) 



P(S>6 C ,5" = 5 C ) 
P(S > 5 C ) ' 
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where 5" represents the linear density on the mass scale M g at redshift z. The relative 
difference between p(X[ = A c , Af> < 0, A3 < 0\5 > 5 C ) and p(5' = 5 C \5 > 5 C ) indicates how 
probable the anti-hierarchical formation of low-mass halos is at given epoch. 

The direct comparison of the two conditional probabilities (eqs.[5] and [8]) is shown 
in Fig.l. The halo mass M g observed at present epoch is set at the dwarf galactic scale 
M g = 1O 6 /i -1 M , and three different cases of the pancake's mass M p are considered: M p = 
1O 8 /i _1 M (dashed); M p = 10 10 h- 1 M Q (solid); M p = 1O 12 /i _1 M (long dashed). As can be 
seen, the value of p{X[ = A c , A' 2 < 0, A 3 < 0\S > 8 C ) is twice higher than that of p(S' = 5 C \5 > 
5 C ) at z ~ 2 for the case of \Q 1Q h~ l M Q <M P < 1O 12 /i _1 M . 

Figure 2 also plots the two probabilities as solid and dashed lines. In this Fig. 2 the 
pancake's mass is set at M p = lO n /i _1 M , and the four different cases of the halo mass 
M g are considered in separate panels: M g = 10 6 /i _1 M Q (upper left); M g = 1O 7 /i _1 M 
(upper right); M g = 1O 8 /i -1 M (lower left); M g = 1O 9 /i _1 M (lower right). As shown, the 
probability distribution p(X[ = A c , A' 2 < 0, A 3 < 0\8 > 5 C ) has a maximum value around 
z = 2, position of which shifts to the low- redshift section as the halo mass M g increases. For 
all four cases of M g at z ~ 2, the value of p(X[ = A c , A' 2 < 0, A 3 < 0\5 > S c ) is consistently 
higher than that of p(S' = 5 C \5 > 5 C ). The results shown in Figs. 1 and 2 imply that the halo 
of mass Mg < lO lo /i~ 1 M observed at present epoch are more likely to have been embedded 
in massive pancakes of mass M p ps 10 u h~ 1 M & around z = 2 rather than formed through 
hierarchical merging. 

Setting the typical mass scale and redshift for the formation of pancakes at 10 u h~ 1 M Q 
and z — 2, respectively, we follow the standard PS approach to evaluate the mass distribution 
function of the low-mass halos that formed anti-hierarchically. According to the theory the 
differential number density of the dark halos in the mass range [M, M + dM] is related to the 
volume fraction F occupied by the proto-halo regions in the linear density field that satisfy 
a specified collapse condition: 
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where A is the normalization factor, which is exactly 2 in the original PS theory (Peacock 
& Heavens 1990; Bond et al. 1991; Jedamzik 1995). If the halos observed at present epoch 
were embedded in massive pancakes at redshift z, the volume fraction F should be written 

as 

F(a) = / d5p(5\X[ = A c , A' 2 < 0, A 3 < 0). (10) 
Js c 

The conditional probability = A c , A' 2 < 0, A' 3 < 0) in this equation (10) can be found 
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from equations (3) and (4) by using the Bayes theorem again: 

P (5, a; = a c , a 2 <o,a 3 <o) 
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Now, the differential volume fraction OF /da in equation (9) is found to be 

/ d5p(5|A c ,A 2 <0,A 3 <0), 
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fd\' 2 f d\' 3 (5 C - A; - A' 2 - A 3 ) p(5, Ac, A' 2 , A' 3 ). (12) 
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The logarithmic derivative of the rms density fluctuation da /din M in equation (9) for the 
case of the sharp k-space filter is also found to be 

<h -^(ln /,-.)• (13) 



dlnM 6a 
where k c is given in equation (1). 

By equations (9)-(13), we evaluate the mass function of the low-mass halos in the mass 
range 10 6 h~ 1 M Q < M < 10 10 /i~ 1 M Q , assuming that all halos in this mass range were 
embedded in pancakes of mass lO n h~ 1 M Q at redshift z — 2. The normalization factor A 
in equation (9) is determined from the constraint that our mass function on the mass scale 
M = lO lo /i _1 M should give the same value as that of the ST formula which is known to 
agree very well with N-body simulation in the mass range M > lO lo /i -1 M0. 

Figure 3 plots our result (solid), and compares it with the original PS (dotted) and the 
ST (dashed) mass functions. As can be seen, our model predicts less number of low-mass 
halos when compared with the PS and the ST mass functions. That is, the formation of low- 
mass halos is suppressed by the earlier formation of massive pancakes. Our mass function 
is found to be well fitted by a power-law, dN/dM w M _L86 , which is shallower than the PS 
and the ST ones, dN/dM xs M~ 2A . This shallow shape of our mass function in the low-mass 
tail is consistent with the recent high- resolution simulation (Yahagi et al. 2004). 
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3. SUMMARY AND DISCUSSION 

By taking into account the possibility that the low-mass CDM halos form through anti- 
hierarchical fragmentation of the massive pancakes, we have derived a new analytic mass 
function for the low-mass halos in the ACDM cosmology with the help of the Zel'dovich 
approximation and the Press-Schechter mass function theory. It has been shown that our 
mass function has a shallower slope in the low-mass tail and predicts maximum five times 
less abundance of dwarf galactic halos of mass 1O 6 /i _1 M than the currently popular Sheth- 
Tormen formula. 

The concept of broken-hierarchy should modify not only the mass function but also the 
other halo statistics from the previous models that were constructed under the assumption 
that the halo formation is perfectly hierarchical. For instance, the two-point correlation 
of dwarf galactic halos would be different in the broken-hierarchy scenario, which in turn 
implies the mass-to-light bias on the dwarf galactic scale would be altered in accordance. 
Our future work will be in the direction of investigating the effect of broken-hierarchy on the 
halo n-point correlations and the mass-to-light bias as well. 

Since our mass function has been derived analytically from first principles without in- 
troducing any fitting parameters, one may not expect it to be very realistic. The formation 
of low-mass halos should be dominated by complicated non-linear processes which cannot 
be described by using first principles alone. However, as it is the first attempt to model the 
broken hierarchy which can accommodate future refinements, it is concluded that our model 
will provide a useful guideline for the theoretical study of the effects of broken hierarchy on 
the structure formation. 

This work is supported by the research grant No. R01-2005-000-10610-0 from the Basic 
Research Program of the Korea Science and Engineering Foundation. 
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Fig. 1. — Probability that a galactic halo observed on mass scale, M g = 10 6 h 1 M Q at present 
epoch was embedded in a pancake at redshift z for the three cases of the pancake's mass: 
M p = Wh- l M Q (dashed); M p = lO^/i^M© (solid); M p = IO^H^Mq (long dashed). For 
comparison, the probability that a galactic halo formed hierarchically at redshift z is also 
plotted (dotted). 
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Fig. 2. — Comparison of the probability that a galactic halo forms in the anti-hierarchically 
(solid) with the probability that it forms in the purely hierarchical way for four different 
cases of the halo mass M g : M g = 1O 6 /i _1 M (upper left); M g = 10 7 /i _1 M o (upper right); 
M g = lO 8 /^ 1 M Q (lower right); M g = 1O 9 /i _1 M (lower left). The pancake's mass is set at 
M p = 10 ll h~ l M & . 
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Fig. 3. — Predictions of our (solid), the Press-Schechter (dashed) and the Sheth-Tormen 
(dotted) models for the number density of dark halos as a function of logarithmic mass. 



